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Abstract. We study the effects of gluon radiation on top production and decay 
processes at an e“'"e“ collider. The matrix elements are computed without any approx¬ 
imations, using spinor techniques. We use a Monte Carlo event generator which takes 
into account the infrared singularity due to soft gluons and differences in kinematics 
associated with radiation in the production versus decay process. The calculation is 
illustrated for several strategies of top mass reconstruction. 


INTRODUCTION 

The study of the top quark is one of the most important goals for future high 
energy experiments. Top has several characteristics which set it apart from the 
other quarks, the most obvious being its large mass: at 175 GeV top is 35 times 
heavier than its partner the b. The large mass of the top quark has a number 
of interesting consequences. It makes ti loop corrections to electroweak processes 
important, so that a good measurement of rut coupled with extremely precise values 
of the W mass will give limits on the Higgs mass. A heavy top may decay into 
supersymmetric (SUSY) particles, providing a way of testing predictions of SUSY 
models. Perhaps most important, large mass means that the top Yukawa coupling 
to the Higgs is large; actually it is close to unity. This means that top studies can 
offer insight into electroweak symmetry breaking, as well as into the fermion mass 
generation process. 

Another consequence of the top quark’s large mass is that its decay width is 
about 1.5 GeV, which means its lifetime is of order 10“^^ seconds. In this short 
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a time, the top quark doesn’t have time to hadronize. As a consequence, the top 
can be completely described by perturbative QCD. Also, the spin information is 
transmitted to the decay products {b and a W^). This allows us to get informa¬ 
tion about the top couplings from angular distributions of the resulting particles; 
differences between measured and SM values would indicate new physics. 

The top was discovered in 1995 at the Fermilab Tevatron proton-antiproton col¬ 
lider. This machine is now being upgraded, and in 2000 the study of the top here 
will start again. Also, at the Large Hadron Collider precise measurements of the top 
parameters (mass, couplings, production cross-section) will be one of the priorities. 

Besides studies at these machines, it would be interesting to study top at a 
high-energy electron-positron collider. The principal advantage of a lepton collider 
versus a hadronic collider is a much cleaner environment; the task of performing 
certain precision measurements will be much easier in the absence of large QCD 
backgrounds. In particular, studies at the tt threshold can be done at a lepton 
collider, but such studies are impossible at a hadron collider. Measurements of top 
quark couplings are especially difficult at a hadron collider. At an e+e“ collider, 
however, besides smaller background, the possibility of using polarized electrons 
in the initial state will be an advantage; also, the coupling of the top to neutral 
currents (7 and Z) will be accessible. In general, the capabilities of a high energy 
lepton collider are complementary to those of a hadron collider. 

To be able to perform these studies, though, a good theoretical understanding 
of top quark’s production and decay processes is essential. An important issue is 
QCD corrections due to real or virtual gluon emission because of the large value of 
the strong coupling constant (~ 0.1). In the following, we will be concerned with 
corrections due to the radiation of a real gluon. 


CALCULATIONAL PROCEDURE 

We are interested in top production and decay with emission of a gluon: 
e+e- ^ 7 *, Z* ^ U{g) hW^W-^g) 

(the top decays into a b — W pair with a branching ratio close to unity). The 
gluon radiation can take place during the tt production process or during the decay 
of either the t or t. Studies of this subject have focused on a single portion of 
the entire process (such as corrections to production [2,3] or decay [4] only), or 
imposed approximations such as soft gluons [5] or intermediate on-shell top quarks 
and massless b quarks [ 6 ]. The results presented here are obtained using an exact 
computation of the entire matrix element for real gluon emission, including all spin 
correlations, top width, and b mass effects [ 1 ]. 

Let us briefly describe the procedure for obtaining these results. We start by 
evaluating the amplitude for the process; there are four diagrams which contribute, 
(see Fig. 1), so we will have: 




FIGURE 1. Feynman diagrams for top quark production with gluon radiation 
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Pt = Pt - ; Pt = Pw+ + Pb (2) 

Ptg = Ptg - ^t + ^'^t^t ; Ptg = Pw+ +Pb+Pg 

and similar definitions for the denominators of t propagators. 

The B, BB terms correspond to the gluon being radiated by the b, b quarks. 
The last two terms correspond to a gluon radiated by the t or f, either in the 
production or decay stage. To separate these cases, we rewrite the propagator 
products as follows: 
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and the amplitude can be writen 
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We identify these three terms as corresponding to gluon radiation in the t decay 
(Al), i decay (A 2 ) or production stage (A3). 

The three separate parts of the amplitude are calculated using helicity amplitudes 
and can be evaluated numerically. We can then compute each of the six resulting 
terms separately: 


^ {|A 1 |^ + |A 2 |^ + |A3|^ + 2 i?e[AlA 2 * + A1A3* + A2A3*]} 

helicities helicities 


To integrate this formula over the phase space we use a three-channel Monte 
Carlo with one channel for each of the diagonal terms |A 1 |^, |A 2 |^, |A3|^ and a 
combination of channels for the interference terms 2 i?e[AlA 2 * -|- A1A3* -|- A2A3*]. 
The phase space region where the gluon energy Eg goes to 0 presents some problems, 
though, because the amplitude has a singularity there. Even with cuts on Eg, the 
rapid variation of the integrand can spoil the integration procedure. To eliminate 
this problem, we tailor the momentum generator to the production of a gluon in 
association with two massive particles ( 7 *, Z* —ttg or t hWg). 


NUMERICAL RESULTS 

In this section, we present some preliminary results obtained using the procedure 
described above. As input, we use the following values: rrit = 175 GeV, T* = 1.5 
GeV, nib = 5 GeV, as = 0.1, and a center of mass energy W = 600 GeV. Also, 
unless otherwise specihed, we impose a cut on gluon energy Eg > 10 GeV. 

First, we note that the radiation of a gluon plays an important role in ti pro¬ 
duction. The lowest-order, tree-level cross section is cto = 0.43 pb, while the cross 
section for the process with emission of a gluon with energy greater than 10 GeV is 
(To = 0.39 pb.^ We see that the two quantities have the same order of magnitude. 

An interesting (and important) issue is top mass reconstruction. This is closely 
related to the problem of assigning the gluon to the correct process in which it was 
radiated. For example, if the gluon was radiated in the production stage, the mass 
of the top is given by 


ml = pI = {pb + Pw+f (4) 

On the other hand, if the gluon was radiated in t decay stage, the above formula 
will give us a low estimate of mp, the correct formula will be 

ml = pig = {pb + Pw+ + Pgf (5) 

Note that this number includes contributions from decay-stage radiation which are not correc¬ 
tions to the total top production cross section. 
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FIGURE 2. Distributions in rht and rhtg as defined in the text. 


In Fig. 2 we present the mass distributions obtained using (4) and (5) respectively. 
While there are clear peaks at the input values of the top mass, we see significant 
tails due to wrong gluon assignments. The existence of such tails can increase 
measurement uncertainties on the top mass or confound attempts to identify top 
events by mass reconstruction. It is apparent that, for good mass reconstruction, 
we need to hnd a method to assign the gluon to the correct process in which it was 
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FIGURE 3. Mass distributions using mass cuts (left) and gluon angle cuts (right). 






radiated. 

Using the variables mt = fhtg = and ifiig = \f^g we define 

four types of events: 

-type 1 : 172 GeU < ifitg^rhi < 178 GeU 
-type 2 : 172 GeU < nit^rhig < 178 GeU 
-type 3 : 172 GeU < fht,mi < 178 GeU 
-type 4 : any other event 

We identify type 1 events as corresponding to gluons radiated in t decay stage. 
In this case we expect pi, ptg to be on shell; however, as we use a Breit-Wigner 
distribution, we accept deviations from the exact value m* = 175 GeV by about 
2ri = 3 GeV. Furthermore, we identify type 2 events as corresponding to gluons 
radiated in t decay stage and type 3 as corresponding to gluons radiated in ti 
production stage. As the gluon energy is quite big (> 10 GeV) compared with F*, 
we expect this interpretation to work; meaning that there will be very few events 
which satisfy one of the conditions 1, 2, 3 but which are open to more than one 
interpretation. Finally, type 4 would correspond to events for which there is no 
compelling evidence for either the production- or decay-stage radiation case; it will 
actually be a mixing of both. 

In the first plot Fig. 3 we present the top mass distribution using events of types 
1, 2, and 3 defined above. We have generated 250,000 events for each case. The 
smooth line is a Breit-Wigner function fitted over the mass distribution. The fit 
reproduces the input parameters with remarkable accuracy; it actually gives us 



FIGURE 4. gluon-b angle distribution 





175.04 GeV for mt and 1.46 GeV for F^. 

Another method for reconstructing the process makes use of the gluon-6 and 
gluon-6 angle distributions. It is known that a gluon radiated by a quark tends 
to go in the same direction as that quark.^ Then, the gluons close to the 6 quark 
probably came from the top decay process; similarly, the gluons close to the 6 are 
likely to have come from t decay. We might guess that the rest of the gluons may 
be assigned to the production process. 

In Fig. 4 we present the distribution of the angle 6hg between the gluon and the 
6 quark. The vertical-hatched part corresponds to type 1 events; and they are close 
to the 6 direction. The horizontal-hatched part corresponds to type 2 events; as 
they will cluster near the 6 quark, and 6 and 6 are mostly in opposite hemispheres, 
these events tend to gather at large angles. The cross-hatched part corresponds to 
type 3 events, and they are distributed uniformly. Finally, the non-hatched part 
corresponds to type 4 events, and, as they are an admixture of the first three, it 
looks like the whole distribution on a smaller scale. 

Using this hgure we can make the following conventions: 

-if Oiyg <0.7 rad assign gluon to t decay 
-if 6j,g < 0.7 rad assign gluon to t decay 
-if 9bg, 6ig > 1 rad assign gluon to tt production. 

With these dehnitions, we construct the top mass distribution presented in the 
second plot in Fig. 3. Again, a £t with a Breit-Wigner function gives us values 
very close to our input parameters: rrit = 175.05 GeV , F* = 1.6 GeV. 

Fig. 5 presents the gluon energy distribution. The soft gluon singularity is 
visible (we have used & Eg > 1 GeV cut). This raises an interesting problem: soft 
gluons are not visible in the detector; therefore, the process in which a soft gluon 
is radiated will give the same signal as the process in which there is no gluon. The 
study of the experimental cross section for this last type of process will have to 
take into account radiation of gluons with energies up to 5 GeV, maybe. This also 
requires including virtual corrections, as in [6]. 

At this point, it is interesting to see if the interference terms contribute. We 
can integrate separately the contribution of diagonal terms and the contribution of 
interference terms; it turns out that this last contribution amounts only to around 
1% of the total distribution, as shown in Fig. 5. 

Up to here, we have worked at the parton level and presumed exact information 
about the momenta of particles in the final state. In a real experiment, though, 
there are a lot of complications: the electrons in the initial state can lose energy 
through radiation before interacting, thus having a lower GM energy; we don’t see 
neutrinos, don’t know which jet corresponds to which quark (or gluon), and the 
energies measured are not exact. One can expect that the detector resolution will 
have one of the biggest effect on the top mass (and width) reconstruction. 

In Fig. 6 we present the distribution for the top mass (using gluon angle cuts) 


Since we keep the b mass, there is strictly speaking no collinear singularity; however gluons 
tend to be radiated from the b at small angles nonetheless. 
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FIGURE 5. Gluon energy distribution; total cross-section and interference contribution, respec¬ 
tively. 



FIGURE 6. Mass reconstruction distribution, taking into account detector resolution (the 
dashed line corresponds to exact energy values). 


obtained after taking into account this effect. The spread in the measured energies 
is parametrized by gaussians with widths a = 0A\/E for quarks and gluon, and 
a = O.lhV^ for the IT’s. The rrit value obtained in this case is still very good 
(175.03 GeV) but the width of the distribution in this case reflects rather the 
detector resolution effects than the top width. 















CONCLUSIONS 


We have presented preliminary results of a calculation of real gluon radiation 
in top production and decay. [1] Mass reconstruction using various mass cuts and 
gluon-bottom quark angle cuts was performed. In an ideal situation, (no initial 
state radiation, perfect particle identification, exact energy and angle measure¬ 
ments) the results were found to be very good, he., the distributions faithfully 
reproduced the input mass and width. In a further study, we will have to take 
into account the effects of energy smearing due to detector resolution, initial state 
radiation, undetected neutrinos and jet identihcation problems. Effects due to 
interference between gluon radiation in the production and decay stages do not 
appear to be experimentally visible, at least for hard gluons. Soft gluon radiation 
(the experimental signal in this case will mimic that of a process without any gluon) 
and virtual gluon corrections have to be studied in more detail. 

We thank C.R. Schmidt and W.J. Stirling for helpful correspondence and dis¬ 
cussions. This work was supported in part by the U.S. Department of Energy and 
the National Science Foundation. 
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